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Abstract. In [f|, we studied host response to a pathogen which uses a cycle of immunologically 
distinct stages to establish and maintain infection. We showed that for generic parameter values, 
the system has a unique biologically meaningful stable fixed point. That paper used a simplified 
model of T-cell activation, making proliferation depend linearly on antigen- T-cell encounters. 
Here we generalize the way in which T-cell proliferation depends on the sizes of the antigenic 
populations. In particular, we allow this response to become saturated at high levels of antigen. 
As a result, we show that this family of generalized models shares the same steady-state behavior 
properties with the simpler model contemplated in Q . 



1. Introduction 

Pathogens that cyclically traverse different stages during their life cycle or during an infection 
process have been studied since the late nineteenth century. Important examples are Plasmodium 
( ID, Tryp anosoma ((l6jp. and the family of herpes viruses, including the Epstein-Barr virus (EBV) 
( Iff, [lj]). One remarkable characteristic of infections with many of such pathogens is life- long 
persistent infection (0,0,[I3],[l,[l6i). 

In Q , we introduced a model of a pathogen that uses a cycle of n antigcnically distinct stages 
to establish and maintain infection. The model is given by 2n differential equations, 

dS 

(1) -± = Fj(S, T) = r , :/, , .V , : - n ,S, - f ,S, - i>,S, l , 

^ = G j (S,T) = c j S j T j -bT j . 

Here Sj denotes the pathogen population at stage j, Tj is the cognate host response. The indices 
j = 0, . . . , n — 1 are taken modulo n. The parameters represent the following processes: 

• aj is the decay rate of stage Sj. If a,j is negative, this state proliferates. 

• fj is the rate at which stage Sj is lost to become (or produce) stage Sj+i. 

• Tj is an amplification factor in the process by which stage Sj becomes (or produces) stage 
Sj+i. For example, the loss of one lytically infected cell may produce Tj = 10 4 free virus. 

• pj represents the efficacy of the immune response Tj in killing infected stage Sj. 

• Cj is the antigenicity of stage Sj . i.e., its efficacy in inducing proliferation of immune response 



b is the natural death rate of the response Tj. We assume it is the same for all stages. 
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We refer to the parameters collectively as 9. Except for aj, j = 0, . . . , n — 1, these are assumed 
non-negative. 

Our flagship result is that while ([T]) has 2™ fixed points for generic values of 9, exactly one of 
these is biologically meaningful and stable 

Let us focus for the moment on the terms —pjSjTj of Fj and CjSjTj of Gj. The term —pjSjTj 
represents the killing of pathogen at stage j (usually infected cells in a particular differentiation 
state) by the cognate T-cell population. This takes place pursuant to an encounter between T-cells 
and infected cells displaying antigen complexed to MHC. To a first approximation the rate of such 
encounters is proportional to the product of the sizes of the two populations. Thus, to a first 
approximation, this term reflects the mechanism of the biological process it represents. 

The term CjSjTj of Gj represents proliferation of T-cells in response to the presence of antigen. 
Here, the underlying biological processes are considerably more complicated, involving a number 
of cell types. Initially, T-cells are activated and begin to proliferate only in response to antigen 
presenting cells, particularly dendritic cells (DCs) (|12|). The density of presented antigen is known 
to affect these T-cell-DC interactions 

(0, EH), eliciting differing CTL responses at different den- 
sities, including T-cell exhaustion at high concentrations of presented antigen ((HI). Activated 
CD8+ T-cells also exhibit central and effector memory phenotypes and the relationships between 
these phenotypes is not well understood ((H). Finally, the length of the cell cycle places a hard limit 
on the rate at which the T-cell population can proliferate. Thus, the rate of T-cell proliferation 
becomes saturated for large amounts of antigen (|9[). In this, they bear a similarity to the rates of 
enzyme catalyzed chemical reactions (reviewed in 2|, see Q for experimental evidence). 
To accommodate dose-dependent effects, we will study the system 



(2) 



dS 

-jj- = Fj(S,T) = rj-xfj-iSj-x - <i,S, - faSj -PjSjTj 
d ^L=G ] {S 1 T) = V3 {S 3 )T 3 -bT r 



which generalizes ([T]). The function H is unchanged from H. The terms of Gj represent proliferation 
of CTLs in response to the presence of antigen and the loss of CTLs due to death or decommis- 
sioning. We use functions ifj : [0, oo) — > [0, oo) to denote the dose-response curves. We assume 
that for each j, <pj(0) = 0, and that for x £ (0, oo), <p'j(x) exists and is positive. In particular, 
each tfj is continuous on [0, oo) and strictly monotone increasing. The possibility of dose-response 
saturation arises from the case where there is an rrij £ M so that lim a; _s. 00 <fj(x) = rrij. We show 
that with appropriate modification, the major results of [f| hold for ([2]). While the term fj(Sj) is 
still phenomcnological in that it omits discussion of biological mechanism, we argue in Section |4] 
that it may well offer a way to address this limitation. 

The system ([2]) can exhibit a behavior which does not arise with (fTJ). In either case if a,j + fj < 0, 
we say that j is self- establishing. It is not hard to see that a self-establishing stage which is not 
regulated will expand without bound. On the other hand, as we will show, if rrij < b, the host 
cannot mount a response to Sj. It is immuno-incompetent with respect to this stage. If j is 
self-establishing and the host is immuno-incompetent with respect to j, we say that the parameter 
set is fatal. If the host is immunologically incompetent at all stages we say the host is totally 
immunologically incompetent. Clearly, in this case, if the basic reproductive number of the pathogen 
is greater than one, infection is also fatal to the host. Accordingly, we will assume that the host is 
immuno-competent for at least one stage. 
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2. Background and definitions 

We start by transforming |2j through a change of coordinates. For this purpose we take rrij := 
ix->oc ^jix) GiU {oo}. Notice that if rrij > b, there is a unique value bj € R so that <Pj(bj) = b. 



If rrij < b, we take bj := oo. So wc define 




if bj < oo 



otherwise 

In the former case, Cj is the marginal antigenicity of Sj at the value bj. We now use the linear 
change of coordinates 

H : R 2n -> K 2 " 
(Sj,Tj) ^ (Sj^j) := := (c^p/Z}). 

This gives the equations 

^ o 

(3) = T) = Tj-xfj-rSj-x - QjSj - fjSj - S^T, 

^=^(S,T)=^)T 3 -&T 3 - 



(4) = Fj(S, T) = Ti-ifr-iSj-t - OjSj - fjSj - S/r, 



Note that for each j, Tpj still enjoys the properties that it is differcntiable, Tp'jix) > for x > and 
<^,-(0) = 0. We now take to be the unique solution to Tpjibj) = b, i.e., bj := Cjbj, where such 
exists. Note that Jpj now enjoys the additional property that Tp'foj) = 1. In the case studied in [f|, 
ifij(Sj) = CjSj and thus Tpj{Sj) — Sj, giving 

dt 

^L=~dj(S,T) = SjTj-bTj 

We will henceforth drop the bars and assume our equations are given in the form ([3]). We take a 
parameter set 6 to be a set of values for b, rj, fj, aj, rrij and bj, j = 0, n — 1. When we need to 
make an explicit comparison with the cpj of ((2]), we will refer to the later as "biological ip" , f] 310 - 

Wc will adopt the following notational conventions. Sets such [j, k) and [j, k) are to be taken 
cyclically. That is to say, if j < k, then [j, k] = {j, . . . , k}, while if j > k, [j, k] = {j, ...,n — 
1,0,..., k}. We take [j,j) to be the empty set so that any product taken over [j,j) is equal to one. 
We abuse notation by taking [0,n) = {0, . . . ,n — 1} 

Wc now review and in some cases generalize the definitions of 

Definition 1. Given a fixed point (S*,T*) of ([3]), the regulated and unregulated stages of (S*,T*) 
are 

Reg(S*,T*) = {j\T;^0} 
Unreg(S*,T*) = {j|T* = 0} 
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Definition 2. (S*,T*) is biologically meaningful if S* > 0, T* > for j = 0, . . . , n — 1. It is 
infected if for some (hence, all, see 1) below) j, S* > 0. 

Definition 3. Given a parameter set 9, the self-establishing stages are 

SE(0) = {j | a 3 + ./' < 0}. 

Definition 4. The immuno-incompetent stages of 9 are 

Incomp((9) = {j \ mj < b} = {j \ bj = oo}. 

Definition 5. If j 6 SE(0) n Incomp(0), we say that the stage j and the parameter set 9 are 
fatal. We will assume that the host is capable of mounting a response to at least one stage, i.e., 
Incomp(6>) ^ [0, n). 

Definition 6. If SE(0) = 0, the follow-on constants of 9 are 

Mj = 

a i+i + Jj+i 

M jk = n M i 

ee[j,k) 

In the case where SE(0) ^ 0, Mj is only meaningful for our purposes for j + 1 ^ SE(0). Accordingly, 
Mjfc is only meaningful if (j, fc] n SE(9) = 0. Note that for every k G [j,£) it holds = Mji-Mki- 

Definition 7. We say that j starves k and write j >- k if bjMjk < bu- Here we assume that Mjfc 
is meaningful and that 6j is finite, though bk need not be. In particular, if Mjk is meaningful, 
j lncomp(0), and k € Incomp(^), then j >- k. 

Definition 8. The starvable stages of 9 are 

Str(0) = {k | there is j so that j >~ k}. 
The unstarvable stages Unstr(0) are the complement of these. 

Definition 9. A biologically meaningful fixed point (S*,T*) is saturate^ if Rcg(5*,T*) = 
Unstr(#). It is moderated if for j € Unreg, S* < bj. 

Definition 10. If SE(0) = 0, we define 

n-l 

R = l[M j . 

3=0 

Rq may be interpreted as the number of copies of the pathogen produced by a single copy entering 
a naive host (j!). It is not hard to see that Rq is invariant under the transformation H as befits a 
property of the organism being described. 

Definition 11. When we say that 9 is generic we will require the following: 

• Ro ± 1. 

• There is no j so that aj + fj =0. 

• There is no pair (j, k) so that bj < oo, bk < oo and bjMjk = bk- 

• At the saturated biologically meaningful fixed point, there are j and k so that T* ^ T£. 

1 There is an unfortunate collision here between the use of the term saturated to denote the host mounting a 
T-response to all stages capable of supporting one and the meaning of the term used in the Introduction above, 
namely, a maximum prolficration rate, with no increase through further stimulation. 
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It is not hard to see that each of these conditions has measure zero, thus justifying the use of 
the term generic. The detailed motivation for these exclusions can be found in 

Definition 12. Suppose Rcg(S'*,T*) ^ 0. Given a stage k, we define hk to be the unique stage 
such that h k 6 Reg(5*,T*) and (h k ,k) C Unrcg(S'*, T*). 

3. Results 

We will start by assuming that SE(0) = 0. This will be a standing assumption until it is lifted 
in Section l3~Tl 

The linear stability analysis performed in [j| is possible because we were able to calculate the 
characteristic polynomial of the Jacobian matrix of the right hand side of ([I]), which corresponds to 
setting ipj(Sj) = Sj for each j in ([3]) (recall that we are omitting the bars). Here we contemplate 
more general cpj : R — > R, j = 0, n — 1 (which is a consequence of contemplating more general 
ifj : R — > R) with the properties mentioned above. Consequently, in order to make use of the 
results obtained in Q, we need to establish what changes are induced on the Jacobian matrix 
through the use of more general functions (fj . The partial derivatives of the right hand side of ( [3]) 
are given by 



r fe _i/ fe _i ifj = fc-l 

-a k - f k - T k if j = k 
otherwise 



dF k _ \-S k if j = fc 
dTj 1 otherwise 

dG k = \y k {Sk)T k itj = k 
dSj |0 otherwise 

dGk = {^k{S k )~b iij = k 
dTj 1 otherwise 

Since the functions Fj do not depend on any ipj, only the partial derivatives 0^ and differ 
from the results obtained in (H- As we shall see, most differences vanish when the functions arc 
evaluated at a fixed point. 

Proposition 13. 

(1) Rq is the basic reproductive number of the pathogen. 

(2) If Rq < 1, the pathogen fails to establish infection and (S*,T*) = is a local attractor. If 
Ro > 1, the pathogen is able to establish infection. In particular, this makes (S*,T*) = 
an unstable fixed point. 

(3) If Ro < 1, (S*,T*) = is a global attractor. 

Proof. There are two ways to establish the first and second claims. One is by using the interpretation 
of Rq in terms of the lifespan and productivity of each stage. The other is by computing the 
eigenvalues the Jacobian. We briefly sketch the first approach. In the absence of immune response, 
stage has an expected lifespan of . ■ During that time, it produces J°+f copies of stage 1. 

These, in turn, produce J°+f B ai+f t C0 Pi cs °f stage 2. Continuing in this way produces Rq copies 
of stage 0. If this is greater than I, the pathogen can establish infection, if less than 1, not. 
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To see these two claims using the eigenvalues of the Jacobian, notice that since <Pj(0) = 0, 
j = 0, n — 1, the Jacobian matrix evaluated at (S*,T*) = (0, 0) is identical to the one obtained 
in [H, Proposition 1]. Thus the claim follows from Propositions 1 and 2 of 

The third claim comes from showing that the reproductive number in the presence of immune 
response is no more than the reproductive number in the naive host as in Proposition 1] . □ □ 

The following correspond to the numbered observations in Section 3 of Q and follow from the 
fixed point equations 

Sj = F,(S* 7 T*) = r, ,./\ ,.s; , - S*( aj + fj +T*)=0 
fj = G 3 {S*,T*) = ((pj(S]) - b)T* = 

1) Given (S*,T*), if there is j such that S* = then (S*,T*) = 0. 

2) If j e Reg(S*, T*), then S* = bj. 

3) If j e Vmeg(S*,T*), then T* = 0. 

4) If j + 1 e Unrcg(S*,T*), then S* +1 = S*M r 

5) If \j + l,k] C Unreg(S'*,T*) then S k = S*M jk . This follows by induction on the previous 
observation. 

6) Assume Reg(S*,T*) ^ 0. If k e Unreg(5*,T*), then 5* = b hk M hhk . This follows from 2) and 
5). 

7) If 6 is generic and (S**,T*) ^ then Reg(S*,T*) ^ 0. Were Reg{S*,T*) = 0, then by 5) 
Sq = SqRq. Consequently Rq = 1, contradicting our first genericity requirement. 

8) If j e Reg(S*, T*) then T* = r -^J^S*_ 1 - {a, + fj). 

9) If j E Reg(S*,T*) then T* = r s ,/, i^-M h , , - (aj + fj). This follows from the fact 
that S*_ x = Sl.Mh^x (which follows from 6) if j - 1 eUnreg(5*, T*) , and holds trivially, if 
j - 1 GReg(5*,T*)) and = b hj . 

10) If j G Reg(S*,T*), then T* > if and only if b hj M hj j > bj. This follows from the previous 
observation; (recall that we have assumed SE(#) = 0). 

Proposition 14. Suppose 9 is generic and (S*, T*) is a biologically meaningful fixed point. Suppose 
further that j >- k . If j £ Reg(S*,T*), then k £ Unreg(S* ,T*). 

Proof. Suppose j>k and j £Rcg(S*, T*). Let {j = j ,j x , ...,j m }= [j, k)(lRcg(S* ,T*) (cyclically 
ordered as listed) and let j m +i = k. If m = 0, then, h k — j and, due to j y k, it holds bh k Mh kk < b k . 
Thus, the claim follows by observation 10), above. Otherwise we have j ^ j m = h k and it suffices 
to show that S* m M Jmk < b k . For I = 0, . . . , m, we have jt eReg(S*,T*) so we must have S* ( = b je . 
For I = 0, . . . , m — 1 we must also have Sj < S* e Mj e j e+1 , because hj e+1 = ji. We then have 

bji = "^j* < S* Mj j 1 
b J m - S* m < S* Q M jojm = b jo M jojm 

so that 

St k M hkk = S* m M jmk < b jo M jok < b k . 
Thus k eUnreg(S*,T*) as required. □ □ 
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Proposition 15. Let 9 be a generic parameter set such that Rq > 1. Then >~ is a strict partial 
order. 

Proof. We must show that >- is anti-reflexive, asymmetric and transitive. The first follows imme- 
diately from the fact that Mjj = 1. 

To see that >- is asymmetric, suppose we have j y k and k y j . We then have bjMjk < b k and 
b k M k j < bj. This gives bj > bjMj k M k j. But Mj k M k j = R , contradicting R > 1. 

To see the third we suppose that j y k and k y I. We consider two cases, k 6 [j, £] and 
£ G [j, k\. In the first case, we have Mj k M k g = Mj#. We then have bjMj k < b k , b k M k £ < bt giving 
bjMjg = bjMj k M k £ < b k M k e < bt as required. In the second case, we have Mj k M k e = RoMj£, so 
that bjRoMjg = bjMj k M k e < b k M k £ < be. Since Rq > 1, this implies j y i as required. □ □ 

Remark 16. Since y is a partial order, it is cycle-free, that is there is no sequence of stages 
jo y jj y ■ ■ ■ y jo- Consequently, we can define the depth of a stage k, d(k) to be the length of the 
longest chain jo y . . . y j^ k ) ~ It follows that Unstr(0) consists of the stages of depth 0. In 
particular, Unstr(0) ^ 0. Note that if 9 is such that no two stages are comparable, then y is empty 
and every stage is ^-maximal, so Unstr(0) = [0, n). Str(#) consists of the stages of positive depth. 
If Incomp((9) ^ 0, lncomp(0) consists of the stages of maximal depth. 

Proposition 17. Suppose that 9 is generic. Furthermore, let (S* ,T*) be a biologically meaningful 
infected fixed point. Then the pathogen populations are moderated at (S*,T*) if and only if the 
immune response is saturated at (S*,T*). 

Proof. We first show that if (S*,T*) is moderated, then (S*,T*) is saturated. 

We claim Unstr (6) CReg(5* , T* ) . Suppose to the contrary j eUnstr(0)nUnreg(S'*, T*). By 
assumption (S*,T*) is moderated, so S* < bj. Since j G Unrcg(S'*, T*), by 7), 6) and 2) above, 
S* = Sfo.Mhjj = bj, Mi, ,. This gives bh^M^j < bj, i.e., hj y j, contradicting the assumption that 
j € Unstr(#). This proves the claim. 

We claim that Str(0) CUnreg(5*, T*). If Str(6>) = 0, this holds trivially. Suppose k £ Str(0). 
Then there is a maximal j so that j y k. Being maximal j € Unstr (9) and thus j €Reg(5*, T*). It 
follows by Proposition [14l that k sUnreg(5*, T*) as required. 

We now show that if (S*,T*) is saturated, then (S*,T*) is moderated. 

If Unreg(S'*, T*) = 0, the claim holds vacuously. Suppose that k SUnreg(S'*, T*). We must 
show that S k < b k . Again, we choose j to be maximal so that j y k. Since (S*,T*) is saturated, 
j EReg(S*,T*) , thus, by 2) S* = b r If [j + l,k) C Unreg(S'*, T*), we are done, for then 
S k = S*M jk = bjMjk < b k . On the other hand if [j + l,k) H Reg(S'*,T*) ^ 0, choose m G 
[j + 1, k) R Reg(S'*,T*) so that m = h k . Since m £ Reg(S* ,T*), by the assumed saturation 
m € Unstr(6 | ). Therefore j )/- to, in other words, bjMj m > b m . If b m M mk > b kl these two 
inequalities would yield bjMjk = bjMj m M mk > b k contradicting j y k. Consequently b m M mk < b k 
must hold. Now we have = S^M mk = b m M mk < b k as required. □ □ 

Theorem 18. Suppose 9 is generic and (S*,T*) is a biologically meaningful infected fixed point 
which is not saturated. Then there is j £ Unreg(S* ,T*) so that for any open neighborhood U of 
(S* , T*) there is a biologically meaningful point x 6 U so that 1 > . In particular (S* , T*) is 
unstable. 

Proof. Since (S*,T*) is not saturated, it is not moderated. Thus, there is j S Unreg(S'*, T*) with 
S* > bj. Since 9 is generic, S* > bj, for otherwise we would have bhjMhjj = bj. In particular, 

bj < oo. It follows that j £ lncomp(0). Since S* > bj, (fj(S*) > b, so -gf L \(S".T*) = Vj(Sj) - b> Q. 
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Let be the unit vector in the Tj direction. Then, for any S > 0, Tj \is*,T*)+Se T . > 0- Thus, in 
any open neighborhood U of (S*,T*), there are biologically meaningful points whose orbits move 
away from (S*,T*). In particular, (S*,T*) is unstable as required. □ □ 

Theorem 19. Suppose that 9 is generic and that (S*,T*) is a biologically meaningful infected 
fixed point. In particular, not all T* are equal. If (S*,T*) is moderated then (S*,T*) is a local 
asymptotically stable equilibrium. In particular, the eigenvalues of the Jacobian matrix J(S*,T*) 
have strictly negative real part. 

Corollary 20. For a generic parameter set, the system (and hence HpJ has a unique biologically 
meaningful stable fixed point. 

Proof. Since the sets of starvable and unstarvable stages depend only on 9, there is exactly one 
saturated fixed point, hence exactly one moderated fixed point. The Corollary now follows from 
the Theorem. □ □ 

Theorem \19\ The proof of the corresponding Theorem in [f| proceeds by showing that the Jacobian 
matrix of the system ((3]) has eigenvalues all of whose real parts are negative. It will therefore suffice 
to show that we can carry out the same computation on the Jacobian matrix of ([3]) evaluated at a 
moderated fixed point (S*,T*). Since H and H are identical, we need only consider the partials of 
G and G. We have 

dG k _ \T h ifj = k 
dSj 1 otherwise 

8G k _ ly k (S k )T k i£j = k 
dSj 1 otherwise 

Now for k 6 Unreg(S*,T*), both of these partial derivatives vanish, while if k £ Reg(S'*,T*), We 
then have S%. = bk so that ip' k (S%) = 1, and once again, the two are identical. 
Moreover, we have 

dGk = is k -b ]£j = k 
dTj 1 otherwise 

dG k _ j(pk(Sk) -b if j = k 
dTj 1 otherwise 

Now if k £ Reg(S'*, T*), we have = b so that S k — b = in the former case, while in the latter case 
we have «S£ = b k so that <Pk(S£) -6 = 0. Finally, in the case where k € Unreg(S'*, T*), the proof 
of @, Theorem 2] appeals to the fact that (S*,T*) is moderated, thus ensuring that Sjl — b < 0. 
Here, the fact that (S*,T*) is moderated implies that < b k so that (p k (S^,) —b<0 and we can 
proceed as before. □ □ 

3.1. Self-establishing stages. We now turn to the case where SE(#) ^ . In this case we need 
the assumption that 9 is not fatal, that is, SE(0) n Incomp(6 l ) = and Incomp(6>) ^ [0, n). 

We start by observing that if SE(0) ^ 0, then the pathogen is viable. Accordingly, in place of 
Proposition [131 we have the following. 

Proposition 21. If SE(9) ^ 0, then (S*,T*) = (0,0) is an unstable equilibrium. In particular, the 
pathogen is able to infect the host. 
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Proof. Suppose that j £ SE(0). Then ^^(o.o) = — a j ~~ fj > 0- This gives orbits with positive and 
increasing Sj inside any open set around (0,0). □ □ 

The numbered observations 1) through 9) listed above hold without change. Observation 10) 
now requires the additional hypothesis that j ^ SE(0), giving 

10') If j 6 Reg(S*,T*) and j <£ SE(0), then T* > if and only if b h .M hjj > b r 

As before, this follows from observation 9). 

Proposition 22. Suppose that SE{6) ^ and (S*,T*) is a biologically meaningful infected fixed 
point. Then SE{9) C Reg(S* , T*). 

Proof. This follows from noting that j e SE(0), S* > and T* = implies Sj > 0. □ □ 
Proposition 23. Suppose SE{6) ^ 0. Then y is a strict partial order. 

Proof. We must show that >- is anti-reflexive, asymmetric and transitive. The first follows imme- 
diately from the fact that Mjj = 1. 

To see that >~ is asymmetric, suppose we have j y k and k y j . This implies (j, k]P\SF,(9) = 
and (fc, j]nSE(0) = 0, contradicting SE(0) ^ 0. 

To see the third we suppose that j y k and k y i. We consider two cases, k 6 \j,£] and 
i € [j, k]. In the first case, we have MjkMu = Mjt. We then have bjMjk < bk , bkMu < bt giving 
bjMjg = bjMjkMki < bkMkt < be as required. The second case would imply (j, k] U (k,£) = [0,n) 
as well as (j, fc]nSE(6») = and (k, £}r\SE{8) = 0, contradicting SE(0) + %. □ □ 

We define Str(0), Unstr(0), saturated and moderated as before. Proposition [TJ] holds in the case 
SE(0) 7^ 0. However, there is a small change in the proof. In the case where SE(0) = 0, we appeal 
to observation 10). In the case where SE(0) ^ 0, we need to note that j y k implies that k £ SE(0) 
and we are thus able to appeal to observation 10'. The proof then proceeds as before. 

The equivalence of moderation and saturation (Proposition I17[) and their necessity for stability 
( Theorem [T5|) can be proved as before, the result of Proposition E^l playing an important role. 

In order to prove sufficiency in the presence of self-establishing stages (Theorem \T§\ and its 
consequences), we rely on Lemma 1 of Section 8 in [B| and the argument provided there after the 
Proof of the Lemma. 

4. Discussion 

In this paper, we have generalized the T-cell activation and proliferation model of Q in order 
to make that model applicable to more realistic antigen dose - T-cell proliferation response curves. 
While in most regimes, we would expect T-cell proliferation to rise in response to increased anti- 
gen, this rate is not driven by the encounter of T-cells with infected target cells, but rather by the 
presentation of antigen to T-cells by dendritic cells ([13). Thus, a model of the mechanisms un- 
derlying T-cell proliferation must include additional cell populations and quite complicated cellular 
processes carried out by those cells ([HI). 

Further, there is a widespread phenomenon that cannot be explained as the fixed point of a 
system like ([3]), namely the existence of long-lived T-cell responses to multiple epitopes, either of 
a single pathogen or in our case to a single pathogen stage. Multiple T-cell responses are often 



2 We take the opportunity to amend Proposition 7 of Section 8 in [f| . The condition Rq > 1 in the statement of 
that proposition is not required, given that we take SE(6>) ^ as a standing assumption for the entire Section 8. 
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modeled as competing for antigen in a predator-prey dynamic. The antigenicity of the T-cell's 
epitope functions as the T-cell's fitness and this leads to a winner-take-all dynamic where the 
response to the most antigenic epitope survives and the others become extinct ( 13 1 , [T3] ) . It seems 
likely that memory T-cells play an important role in the survival of multiple responses. However, 
the interactions between effector and memory populations are still not well understood ([3], (H). 
These populations exhibit differing longevity. Stated in terms of our model, they do not share 
a common value for b. In addition, high levels of antigen can lead to CTL exhaustion ([ll[), a 
phenomenon that argues against a monotone increasing tp. 

The upshot of this is that if we wish to refine the cyclic pathogen model to present an increas- 
ingly detailed picture of the cell populations and their mechanisms, we will need to include multiple 
immune cell populations at each stage. The dynamics of such a system could be quite compli- 
cated. However, there is a variable that summarizes the collective immune pressure against a given 
stage, namely their net kill rate of the effector populations. Thus, if Tj\, . . . ,Tjk are the effector 
populations against stage Sj, we can write Tj = Tj(Tji, . . . , Tjj~) so that we now have 

= r j ifj :-S i - ( a j + Si + T -j) S r 

We cannot expect that the dynamics of such an expanded system can be mapped to the dynamics 
of ([3]) because we cannot necessarily expect Tj\, . . . , Tjk (and any non-effector populations) to vary 
in a way which makes a function of Sj and Tj. However, once the dynamics of these populations 
are understood, in the neighborhood of a fixed point, understanding the marginal response of Tj to 
Sj may allow us to use to summarize these dynamics in a way which will allow us to establish 
the existence of a stable fixed point. 
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